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ABSTRACT 

Following the recent suggestion that it is dwarf galaxies in clusters - as opposed to large 
^ . ellipticals - that provide the intracluster gas, we estimate the metallicity of the intracluster 
medium (ICM) in such a case. We derive analytical expressions for the fraction of mass of 
Q\ ■ dwarf galaxies that is ejected, and estimate the metallicity of the resulting intracluster gas. 

We find that the metallicity resulting from this hypothesis is adequate only for clusters 
with low-metallicity gas. Since galactic winds from dwarf galaxies occur much earlier than 
those from ellipticals, due to the smaller binding energy of the gas, we show that the gas 
ejected by dwarf galaxies is enriched mostly by Type II supernovae, as opposed to Type I 
j> ! supernovae in the case of large galaxies. We further point out that the gas in small scale 
^ . structures, which never cools and forms stars due to low temperatures and, consequently, 
large cooling time scale, gets dispersed in the process of hierarchical clustering, and is 
incorporated as the diffuse intracluster gas. We show that this process can provide enough 
hot gas to be compatible with X-ray observations in rich clusters. 

Subject headings: Galaxies : abundances - Galaxies : clustering - Galaxies : evolution - 
Galaxies : intergalactic medium 



1. INTRODUCTION 

The existence of hot gas in clusters of galaxies was established from X-ray observations 
almost two decades ago, and so was its metallicity, especially from the emission lines of 
Fe, Si, S, Mg and O (e.g., Sarazin 1988). The amount of the hot gas has been inferred to 



be a substantial fraction of the total gravitational mass of clusters (e.g., White et al.1993). 
For the origin of this gas, several hypothesis have been put forward, the most pursued 
one being that of galactic winds from elliptical galaxies driven by supernovae (Arimoto & 
Yoshii 1987; Matteucci & Tornambe 1987; David et al.1990). 

Several authors have, however, found that the amount of gas expelled from the vis- 
ible ellipticals and SO galaxies is not adequate to account for the total gas mass (David 
et al.1990; Matteucci and Vettolani 1988; Okazaki et al.1993). A large portion of the ICM 
gas, even as large as 90% of it, has therefore been thought to be primordial. This inter- 
pretation, coupled with the fact that the ratio of the intracluster gas mass to the stellar 
mass in the galaxies is higher in rich clusters, has led to the idea that galaxy formation is 
less efficient in rich clusters (David et al.1990; Matteucci and Vettolani 1988). 

The metallicity of the ICM gas has been estimated by various authors from the hypoth- 
esis that the metals are produced and expelled in galactic winds from ellipticals. Chemical 
evolutionary models of ellipticals have been used to predict the metallicity (e.g., Arimoto 
& Yoshii 1987; Matteucci & Tornambe 1987). Okazaki et al.(1993) found that these models 
are adequate to explain the observed total mass of iron in the Virgo cluster, but not to 
explain the total gas mass. Arnaud et al.(1992) defined a dimensionless ratio of iron mass 
to that of the stars in the galaxy, to analyze the data, and came to the conclusion that one 
needed a bimodal star formation in cluster galaxies. Renzini et al.(1993) found this idea 
rather ad hoc, and, introducing another ratio, that of the iron mass to blue light, they 
found constraints from the data on the yields of Type I and II supernovae rate that may 
explain the observations. It seems, therefore, that no single model satisfactorily explains 
the iron mass (or the metallicity) and the total gas mass of ICM, simultaneously. 

Recently, Trentham (1994) has suggested that if dwarf galaxies, at the lower end of the 
luminosity function, are numerous, i.e., if the index of the Schechter luminosity function 
a £j — 1.7, as is seen in some clusters, then they may supply enough gas. These galaxies, 
with shallow gravitational potential wells, and with small binding energy of the gas, are 
vulnerable to huge mass losses through winds driven by supernovae explosions from the 
first generation of star formation (Larson 1974; Saito 1979). Dekel & Silk (1986) showed 
that galaxies with velocity dispersion less than ~ 100 km s _1 would lose much of their gas 
due to large galactic winds, and become low surface brightness dwarf galaxies. Trentham 
(1994) argued that if a fraction 7 ~ 0.05 — 0.33 of the total mass of the galaxy is expelled, 
and if the index of the Schechter luminosity function a ~ ( — 1.9) — ( — 1.7), for masses 
10 4 — 10 11 M , then the amount of gas mass expelled can account for the total gas mass 
in the ICM. 

One of the tests of this interesting model lies in the resulting metallicity of the ICM 
gas. Since the galactic winds from dwarf galaxies (of total mass ^ 10 11 M ) occur in a time 
scale t ^ 10 9 yr, Type I supernovae can not contribute substantially to the metallicity. We 
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suggest here that the hypothesis of dwarf galaxies as sources of ICM gas leads necessarily to 
the dominance of Type II supernovae as the origin of the metals in ICM. In this paper, we 
estimate the total metallicity of the ICM resulting from galactic winds from dwarf galaxies 
using more physical values of 7, rather than a universal number for all galaxies, and we 
derive the relations between the metallicity and various parameters of star formation. 

The structure of the paper is as follows: In §2 we calculate the galactic wind epoch for 
dwarf galaxies and the mass of iron that is contributed by individual galaxies. In §3, the 
total mass of iron in the intracluster gas is estimated using the mass function of galaxies. 
We discuss the nature of the metallicity of ICM in §4. In the paper, we write the present 
day Hubble constant as Hq = 5O/150 km s _1 Mpc -1 . 

2. GALACTIC WIND FROM DWARF GALAXIES 

Galactic winds are thought to be initiated when the thermal energy of the gas in the 
parent galaxy exceeds its binding energy. Below, we will consider galaxies with masses 
below Mr ~ 10 11 M , where Mr denotes the total mass. We will assume that a fraction 
(1 — g) of the mass is in the form of dark matter, i.e., the mass of gas (M g ) and that of 
stars (M a ) comprise a fraction g of the total mass. A few million years after the onset of 
star formation, the first generation of supernovae begins to explode (the life time of a 100 
M star is ~ 3.2 x 10 6 yr), and they put energy into the interstellar medium. 

The energy in the supernovae blastwaves is turned into the thermal energy of the 
gas, as the supernovae remnant (SNR) shells overlap and the thermal gas inside the shells 
interacts with the ambient gas in the ISM. (The energy in the dense shell is quickly radiated 
away and does not contribute much in pumping energy into the ISM gas.) The radii of 
SNR expand as (Cox 1972) 

RsNR,(t) = Ro(t/trad) 2 ^ , for t < t rad 

= Ro(t/trad) 2/7 , for t>t rad (1) 

Here #0 = 0.039 n~ 0A kpc and n is the particle density. The total energy liberated in 
a typical supernovae is assumed to be 10 51 ergs. Half of the thermal energy of the SNR 
is radiated away after a time t raa \, which, for primordial gas, is approximately equal to 
1.7 x 10 5 n -0 - 49 yr for a primordial gas (Babul & Rees 1992). For a mildly enriched gas, 
trad ~ 0.9 x 10 5 A~ 5/17 n~ 9/17 (using a cooling rate of 3 x 10~ 18 AT _1 n 2 ergs cm -3 s -1 , as 
in Dekel & Silk 1986), where A represents the metallicity of the gas (A = 1 for Z = 0.01, 
and higher for larger value of Z). Since these two time scales are within a factor of order 
unity, we will use the one for a primordial gas for simplicity. 
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If we define (as in Babul & Rees 1992) t snr , the time scale for the overlapping of 
SNRs, as the time when the remnant shells occupy ~ 60% of the total volume, then t snr 
can be estimated in the following way. We use for the radius of the galaxy, the expression 
derived by Saito (1979) (or. e.g., Matteucci & Tornambe 1987), which is valid over a large 
range of masses, from globular clusters to large ellipticals, viz., 

= / M T \ °- 55 

where Mt is the total mass of the elliptical galaxy. (This is close to the estimate of the 
radius of the starburst region of dwarf galaxies in Babul & Rees 1992.) We assume a 
constant star formation rate for analytic simplicity, and write the star formation rate as 
M* ~ (rgMx/t / /), where tff = 1.9 x lO 7 ^ 1 / 2 ?? -1 / 2 yr is the free fall time of the gas, and 
r(< 1) is the efficiency of star formation. We will see later that the final results of the 
metallicity of ICM do not depend on r. The rate of supernovae is given by the v = 10 _35 z/5o 
per gram, where z/50 refers to one supernova per 50 M (see also, Dekel and Silk 1986). 
Therefore, the number of supernovae in time t is given simply by NsN(t) = vM*t. 

The rate of supernovae, as denoted by v, of course depends on the slope of the initial 
mass function (IMF). For an IMF with slope x [4>{m) cx m~ x dm\, the factor z/ 50 can be 
written as, 

where M/ and M u are the lower and upper limits of masses of stars in the main sequence, 
and Mi jSn is the lower limit of the mass for the progenitor star for a Type II supernova 
(all in the units of M ). For Mi = 0.1, M u = 100, and M t}Sn = 8, we get u 50 ~ 0.37 (for 
x = 1.35, Salpeter IMF), and z/ 50 ~ 0.90 (for x = 0.95). The term z/ 50 , therefore, succinctly 
expresses four important parameters (x, M;, M u , M^ sn ) and it will be useful in comparing 
our results later with the numerical results of previous workers. 

The timescale of the SNR collision, t snr , can then be estimated as (in the units of 10 9 

yr), 

/ 7/ C n \ -0.54 / n \ -0.27 / \ 0.35 

t snr ^6xlO- 3 n°- 13 (^M ( f T ) Gyr. (4) 

V0.35/ V0.1/ V1O 12 M / y K ' 

Larson (1974) showed that the fraction of the energy of the blast (= E$n) that is finally 
imparted to the gas is 1] ~ 0.1 (see also Babul & Rees 1992), when the remnant shells 
overlap. We will see below that t snr in eqn (4) is much shorter than the time scale over 
which the galactic winds are excited. Therefore, the assumption that 77 ~ 0.1 is valid, since 
the whole galaxy would be pervaded by SNRs by the epoch of galactic winds. 
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As the continuing SNR collisions heat up the gas, the thermal energy content of the 
gas (E t h) increases as, 

E th (t) = 0, t<3.2xl0 6 yr, 

« 7]E SN vM*t , i»3.2x 10 6 yr . (5) 

The epoch of the onset of galactic wind, t gw , is defined as the time when the thermal 
energy of the gas exceeds the binding energy, i.e., Eth(t gw ) > E(,i n {tgw)- [We will neglect 
the effects of stellar winds here (see Gibson 1994).] We take the binding energy of the gas 
as (Saito 1979; also, eqns (1) and (18) of Arimoto and Yoshii 1987 for g <C 2), 

As noted by previous workers, the energy accumulated in the gas in dwarf galaxies exceeds 
the binding energy of the gas even when a small fraction of the gas is turned into stars 
(M a < M g ) (Larson 1974; Saito 1979; Dekel & Silk 1986). We have accordingly assumed 
here that M s £ M g . Eqns (5) and (6) then give, 

V -8.2 x 10- 3 n-VV/ 2 f^)"Y^)" 1 (^)- 1 f^^)°- 45 Gyr. (7) 
9 y V0.1/ V0.35/ V 0.5 ; V1O 9 M / y v ; 

This time scale can also be written as 

V - 9 x ID" 3 ( JL) - 1 (2*°_)-i (Z_) - 1 (Wr_) °™ G yr , (7a) 
9 V0.1/ V0.35 ; V0.5/ V1O 9 M / y ' v ; 

using the density of gas as indicated by the radius from eqn (1). Note that t gw given 
above is the time taken by SNRs to excite the galactic wind, after the first generation 
supernovae explode, i.e., after ~ 3.2 x 10 6 yr. Therefore, to compare with the values 
obtained numerically by previous workers, e.g., Yoshii & Arimoto (1987, Fig.l), Matteucci 
& Tornambe (1987, Table 1), one should add ~ 3.2 x 10 6 yr to the t gw from eqn (7), to 
find the age of the galaxy at the time of the galactic wind. The numbers from eqn (9) and 
those in the previous numerical works are consistent within a factor of a few. However, 
eqn (7) should not be taken too literally for larger masses (M T >, 10 11 M ), for which the 
assumption of constant supernova rate fails. But for dwarf galaxies (M T < 10 11 M ) this 
assumption is a reasonable one, since the time scale of the galactic wind is much shorter 
for these galaxies than the time scale of the decline of the supernovae rate after the first 
peak (Matteucci & Greggio 1986). 

It can also be seen from the supernovae rates of Type I and Type II, as computed by 
Matteucci & Greggio (1986), that SN Is do not begin to contribute substantially to the 
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enrichment of the gas until after a few times 10 8 yr. This is larger than the galactic wind 
time scale for galaxies of our interest. Therefore, to calculate the metallicity of the gas 
expelled by the dwarf galaxies in clusters, we will only consider Type II supernovae. 

Renzini et al.(1993) compared the amount of iron produced in SNII, as predicted by 
various models, and found that an average amount of ~ 0.07 Mq of iron is produced per 
SN II, and that it is fairly independent of the slope of IMF. We write the total amount of 
iron in the galaxy at the time of the galactic wind as, Mp e ~ 0.07z/M*t 3TO , using t gw from 
eqn (9). Approximately, the stars that form within ~ t gw after the onset of star formation 
are the ones that become supernovae before the gas is ejected out of the galaxy, and only 
those stars get to enrich the gas in the galactic wind. This yields, 

M Fe ~ 0.07uMJ gw (M Q ) 

V0.1/ VO.1/V1O 9 M / ° v ; 

The average metallicity of the galaxy at the epoch of galactic wind, is then given by 
[dividing eqn (8) by gMr, since the total content of iron would be shared by gas and the 
stars] , 

0.056 



p-rY-^Y 0) 

\0.l) \W 9 M Q J ' 1 } 



We have used Z Qj p e = 1.89 x 10~ 3 (Renzini et al.1993). The values predicted by the above 
equation match well those from Fig.l of Yoshii & Arimoto (1987) for galaxies M T £ 10 11 
M (see also Arimoto & Yoshii 1987, Matteucci & Tornambe 1987). However, one should 
remember the different assumptions adopted by these authors when comparing the results 
(e.g., the limits on the masses of supernovae progenitors, slope of IMF etc.). We find that 
the values given by our equations are within a factor of a few from previous numerical 
results. Moreover, the equations above are more transparent with all the dependences 
shown explicitly, than the numerical results. 

We can estimate the gas mass in the wind following Larson (1974). We define the 
gas mass and mass which has formed into stars, as M g and M s , respectively. The binding 

E ( \0.45 

energy of the gas is 9 ^ in ~ 1.5 x 10 47 ( 10 1 ^ Mq J erg Mq 1 [see eqn (6)] . The energy in the 
gas can be estimated as 10 48 (^ 5 |)(^-)M S erg, for a typical supernova energy of 10 51 erg. 

M ( \-0.45 

This means that gas will escape from the galaxy when jf- ~ 6.67 f 10 a^ J (^ s f)((n") ~ 
( m \-o.45 

( 6 78x101° m q ) (fi^Xo 2 !)- Therefore, remembering that gMr = M g + M s , one can 
write the fraction of mass that is expelled from dwarf galaxies, as 



, = ^. fi i (_?_\ U.78xlQi°M J WW ( . 

mi xu.j./ i -|- l 6 . 78xl0 io M ) lo.5Jvo.lJ 
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However, this expression for 7 is incorrect at the low mass end, because it predicts 
an unrealistically large value. The above expression neglects the fact that a lot of mass is 
locked up in low mass stars at the time of the explosion of the first generation supernovae, 
and it underestimates the mass of stars at the time of the galactic wind. To take this 
into account, we first note that the galactic wind for galaxies smaller than gMr ~ 10 7 
M© occurs soon after the first generation of supernovae explode. We then note that the 
fraction, mi, of the stellar mass (M s ) that is in low mass stars, can be written in term of 
the IMF slope x, as 

/ 1 1 \ / 1 1 \-i 

For Mi, M u , and M Z)Sn = 0.1, 100, and 8 M , respectively, mi = 0.86 (x = 1.35) and 
mi = 0.59 (x = 0.95). Now, the gas mass in the galaxy decreases due to star formation as, 

where we have used the earlier expression for the star formation rate, and used a value of 
n that is given by eqn (2) for ionized gas. The mass of gas that remains in the galaxy after 
t ~ 3.2 x 10 6 yr, the epoch of the first generation supernovae (~ life time of very massive 
stars) is given by, 

( ( M \-°- 325 \ 

M;= 9 M T exp(-0.16rUwy )' (13) 

However, this neglects the fact that by the time of the first generation supernovae, the 
massive stars [of total mass, (1 — mi)M s ] would return most of their mass back into the 
galaxy. Therefore, the actual gas mass M g = M' g + (1 — mi)(gMT — M' g ). This gives us an 
expression for 7 that is more realistic for low mass galaxies, 



1 = 9 



I - mi 



1 — exp — 0.16r 



/ Mr 
\10 9 M Q 



(14) 



We will use a value of 7 that is smaller of the two values predicted by eqns (10) and (14). 

The resulting 7 is plotted against galactic masses at the present time, M = My(l — 7), 
in Fig. 1(a). It shows a maximum around M ~ 10 8 ~ 9 M & : a large amount of gas is lost 
from such typical dwarf galaxies. For larger galaxies, the deeper gravitational potential 
does not allow the hot gas to escape, while for smaller galaxies, most of the gas is converted 
into stars prior to galactic wind. These values of 7 are physically more realistic than the 
constant 7 for galaxies of all masses used by Trentham (1994). 
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Figure 1(b) shows the corresponding iron abundance of expelled gas in units of the 
solar value. Smaller galaxies with M ~ 10 10 M© have less time to enrich the gas (t gw ^ 10 8 
yr), giving a low metallicity, of order less than 0.5 Z & . If the conversion efficiency of 
the blast energy into the gas, denoted as rj, is lowered, the epoch of galactic wind is 
prolonged and thus the metallicity is increased. The mass of expelled iron is estimated 
as 'jMtZfc = Y^MZFe, and the resulting iron masses for M ~ 10 9-11 Mq are in good 
agreement with those numerically derived by Matteucci & Tornambe (1987) and David 
et al.(1990). 



3. HOT GAS AND IRON IN COMA CLUSTER 



Having derived the physically plausible expressions for the ejected gas fraction 7 
of initial total mass and its iron metallicity Zp e in each dSph galaxy, we estimate 
the total amount of hot gas and iron in the ICM. Following Trentham (1994), we use 
a Schechter luminosity function (j>(L) with the faint-end slope of a, i.e., 4>(L)dL = 
4>* exp(^— j^j (jfr) 7jt- Also, as in Trentham (1994), we use a scaling law for the mass- 
to-light ratio, (M/L) oc L^, for dSph galaxies. The characteristic luminosity L*, at the 
knee of Schechter function, is ~ 1.42 x 10 9 /i5 2 L© for dSph galaxies (Sandage et al. 1985). 
Converting into the mass function with M* /L* = 15 for the characteristic mass (Trentham 
1994) and using 7, we derive the total mass of gas ejected from dSph galaxies as 

irM/ » f M+ 7 / / M \V(/3+i)\ / M dM 
M gas = N M / ^-exp - — — — —^ttj- , (15) 



M_ 



7 \M*J J\M*J (l + p)M 



where M denotes the present-day galaxy mass, and A^* is a characteristic number of 
dSph galaxies, which is equal to 302.6 for Coma (Trentham, private communication). The 
integral is performed over all masses of dSph galaxies between M_ and M + . Similarly, the 
final metallicity of the cluster gas is estimated as 



F "~ !^(M)dM ' <16) 

where Z Fe is estimated from eq.(9). 

Adopting the mass of observed hot gas in Coma within 5h^Q Mpc as ~ 5.1 x 
10 14 /i^ 2 - 5 M G (Briel et al. 1992), we derived the required combinations of (a, j3) to re- 
produce all the gas in Coma by galactic winds, and plotted them in Fig. 2(a). The thick 
solid line corresponds to a low-mass bound of M_ = 1O 4 M , while the thin line, to 
M_ = 1O 6 M . We adopt h 50 = 1 here, but different values, such as h 50 = 2 give es- 
sentially the same results. Comparing with the estimate by Trentham (1994) with the 
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assumption of 7 = const., our physically realistic values of 7 imply more stringent bounds 
on (a,/3); when M_ = 1O 4 M and x = 0.95, our required combinations of (a, (3) are sim- 
ilar to his minimum case of 7 ~ 0.065 = const., and a Salpeter IMF and/or larger M_ 
give smaller values for (a,f3). For instance, if the observationally allowed lower limit for 
the faint-end slope of the Schechter function is —1.9 < a < — 1.7 (Driver et al. 1994), 
Fig.2(a) implies -0.55 < (3 < -0.37, compared to -0.42 < (3 < -0.27 in Trentham 
(1994). Thus it becomes more difficult to reconcile with the observationally derived value 
of (3 = -0.22 ± 0.09 for dSph galaxies (Kormendy 1990). 

Figure 2(b) shows the resulting iron metallicity in Coma, in the hypothesis of dwarf 
galaxies as sources of the intracluster gas. The thick solid curve presents the case with 
parameters (M_ = 10 Mq,tj = 0.1), indicating the very low abundance below 0.01 com- 
pared to the typically observed value ~ 0.2 (Hughes et al. 1988). Changing M_ from 
1O 4 M to 1O 6 M (so that the gas ejection is prefrentially from more massive galaxies) 
or from 77 = 0.1 to 0.05 (so that t gw is prolonged), simply doubles the iron abundance as 
0.01 ~ 0.025 and does not produce the observed high abundance. Thus we conclude that 
even if the number of dSph galaxies can supply the whole cluster gas, they are not the 
main sources of high metallicity observed in many clusters, because of the rapid loss of gas 
after the first burst of star formation. The smaller the low-mass limit of dSph galaxies, 
the smaller is the iron abundance, due to the increased contribution of smaller galaxies to 
ICM. The effects of Type la SNe, even if they contribute, do not provide more than three 
times iron of the value obtained here (Renzini et al. 1993). Therefore, we conclude that 
galactic winds from normal bright ellipticals are the main sources of intracluster metals 
(e.g., Matteucci & Vettolani 1988; Arnaud et al. 1992; Okazaki et al. 1993). 

4. DISCUSSION 

Here we discuss some implications of our results above. First, we discuss the origin of 
the intracluster gas, in the light of the constraints on the luminosity function of dSphs in 
clusters to be able to supply the hot gas, as plotted in Fig. 2(a). We consider sources for 
the primordial gas that can explain the X-ray observations. Then, in §4.2, we discuss the 
implications of our results on the iron abundance of the ICM in the case of dwarf galaxies 
being the dominant contributors. 

4.1. Origin of the ICM 

It has been argued by some previous workers that most of the X-ray emitting hot gas 
in clusters ought to be primordial, because the total gas ejected from bright ellipticals is 
insufficient to account for the large mass of hot gas. Okazaki et al.(1993), for exmaple, 
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estimated that less than 10% of the observed gas mass can be explained by gas ejected from 
large ellipticals. The primordial gas may have pervaded the cluster from the beginning, 
supplied by later gas infall (White et al. 1993), or gas which failed to accrete into local 
galactic potentials (Evrard et al. 1994). In contrast, Trentham's (1994) idea is that dwarf 
galaxies can supply most of the cluster gas via galactic winds, on the assumption that 
the fraction 7 of ejected gas mass is constant for all galactic masses. Here, based on more 
realistic estimates of 7, we have found that in Coma, one needs smaller (a, (3) than his case, 
i.e., more unlikely values than allowed in previous observations (a ^ — 1.8,(3 ~ —0.22). 
Although measurements of exact values of (a, (3) still elude us, this hypothesis seems 
difficult to support. Moreover, as we shall explain below, gas in galaxies with M < 1O 6 M 
may not cool to form galaxies, and this makes the situation worse (Fig. 2a). 

In the context of CDM theory, most of the baryonic gas is condensed into a number 
of small objects; as a matter of fact, the number can be very large due to the steep slope 
of mass function ~ —2. However, the temperature of small-scale virialized gas is below 
10 4 K for all M ~ 1O 6 M and < la density contrasts of M ~ 1O 7 " 8 M , for which the 
cooling rate of the gas via radiation is essentially null (in the absence of molecular gas). 
Lacey & Silk (1991) showed that protogalactic objects with mass below M ~ 1O 6_7 M 
never cool within the Hubble time, and this lower limit of the cooled mass is not changed 
in the presence of photoionization by diffuse UV background radiation (Chiba & Nath 
1994). Probably such very small and dark objects are observed as Lyman a clouds in the 
intergalactic space, and the gas is stably confined in mini CDM halos (Rees 1986), with 
the temperature being kept ~ 10 4 K in the presence of photoionization (Efstathiou 1992). 
Therefore, there may be a number of small objects whose baryonic components have not 
cooled and that thus have failed to form stars. 

We remark here that according to numerical simulations of cluster formation (see 
the review by White 1995), small-scale structures, which were present in earlier times, 
disappear due to mergings and tidal forces in the progress of hierarchical clustering, and 
the final output is a regularly structured cluster with a smooth profile. It is thus likely 
that uncooled gas which was originally condensed into small objects is dispersed during 
hierarchical clustering, and heated up to the virial temperature of the cluster after the 
collapse of the cluster as a whole. 

If this is the case, we should add this uncooled primordial gas which is originally 
resident in small scale structures, M uncoo i ed , to the cluster gas. We estimate this as, 



where Mn m denotes the lower limit for the mass of virialized objects. The dotted line in 
Fig. 2 (a) corresponds to the case of Mu m = 1O 2 M with M_ = 1O 6 M , and is comparable 




(17) 
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to the case of Trentham's upper limit, 7 = 0.33 = const, (the maximum baryon fraction 
in Coma). We note that this uncooled gas can account for about 90% of total cluster gas 
in this case, and more than 80% even when Mi im = 1O 4 M as plotted in Fig. 3. Therefore, 
in order to reconcile the idea of dSphs supplying the intracluster gas, with the present 
observational lower limits for (ct,/3) based on realistic estimates of 7, it is important to 
consider the non-negligible fraction of uncooled primordial gas in the mass scale below 
~ 1O 6 M . Furthermore, we note that the slope of the corresponding galaxy mass function 
(j)(M), given as (a— is ~ —2.2 for Mn m = 1O 2 M , which is roughly in agreement 
with CDM slope of ~ — 2. On the other hand, without uncooled primordial gas, the slope 
is ~ -2.4 for M_ = 10 4 M o and ~ -2.6 for M_ = 1O 6 M , which seem to be too steep. 

Therefore, we propose the uncooled gas in small objects as an alternative candidate 
for the source of the primordial gas that can explain most of the intracluster gas. This may 
also agree with the observed trend that richer clusters hold more hot gas; small objects 
below 1O 6 M may be more vulnerable to destruction in richer clusters, thereby providing 
more gas. It is however premature now to test this hypothesis from numerical simulations, 
because the existent calculations lack enough dynamic range (e.g., a gas particle mass in 
Evrard et al. 1994 is ~ 1O 8 M ). More sophisticated simulations including better resolution 
are required to check this possibility. 

We note that Trentham (1994) used a universal value of a throughout the cluster, 
which need not be true. The dwarf galaxies near the central region of the cluster are more 
vulnerable to being torn apart by tidal forces than the dwarfs in the outer region, and 
this will mean that a varies with radius in clusters. It is possible that the dwarf galaxies 
contribute to the ICM gas more in the outer regions than near the central region. Babul 
& Rees (1992) have argued that the inhibition of galactic wind from dwarf galaxies in 
clusters ( after the ICM has formed) may give rise to some observed differences between 
dwarf galaxies in the central regions and those in the outer regions of a cluster. 

4.2. Metallicity of ICM 

In the previous sections, we have argued that the metallicity of the ICM gas originates 
mostly in Type II supernovae, if the dwarf galaxies are indeed the sources of the observed 
hot gas. We suggest that the ratios of different elements of the ICM gas should provide 
a test of the hypothesis of the dwarfs as main contributors of the gas. It is interesting to 
note that some previous workers (e.g., Canizares 1988; White 1991) have already inferred 
the origin of the metallicity of ICM in Type II supernovae, from the observations of the 
ratio of iron to oxygen. The observed ratio of [O/Fe] is known to be (3 — 5) times the 
solar value, which is indicative of Type II supernovae. This is in contrast to the prevalent 
notion, from the models of galactic winds from large ellipticals, that the bulk of the iron is 
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produced in Type I supernovae, where the oxygen to iron ratio is much smaller than that 
in Type II supernovae. For example, the ratio [O/Fe] predicted by Matteucci & Vettolani 
(1988) is as low as ~ 0.9, and in other models (e.g., David et al.1990) it is of the order 
~ 1.5 (Okazaki et al.1990). 

The dominance of Type II supernovae as the contributors to the metallicity of the ICM 
gas can be increased by further interactions between the galactic winds and the intracluster 
gas. It is possible that, if the explosions that ejected the gas out of the galaxies happen 
at the same epoch, then the ejected gas, now in the ICM, may inhibit further galactic 
winds from any dwarf galaxies that may have retained some fraction of its gas to undergo 
starbursts later. The heating time scale of the intracluster gas depends on the mechanism 
of heating (Sarazin 1988). If the ejected gas is heated in a time scale shorter than that 
of the later starbursts from the surviving dwarf galaxies, then, as Babul & Rees (1992) 
argued, the winds can be confined by the surrounding hot gas. They estimated that, if 
(nT)icM ~ 10 4 cm -3 K, then the wind from galaxies of mass My ^3x 10 9 M would 
not extend beyond the galaxy's dark halo. The gas in the wind in this case would cool and 
probably fall back on to the core. It is possible that this effect may inhibit the injection 
of metals from Type I supernovae, which begin to enrich substantially the parent dwarf 
galaxies after a few times 10 s yr of the first generation of stars. 

Since our estimates of the mass of iron of the ICM from dwarf galaxies fall short for 
clusters with large metallicities, it is necessary to take the large ellipticals into account in 
some cases. Although the importance of Type I supernovae in enriching the gas inside and 
in the galactic winds from large ellipticals has been discussed in the literature, it is still a 
debatable issue. It is possible that the gas inside massive galaxies and the gas in galactic 
winds are enriched by different sources. Detections of any signature of Type II supernovae 
origin of the ICM gas will not necessarily rule out contributions from large ellipticals, if 
Type I supernovae are found not to be important for enrichment in the galactc wind gas. 

It is possible that in rich clusters, where tidal forces are stronger, the dwarf galaxies 
do not supply most of the gas in ICM, and the metallicity is larger because the intracluster 
gas is mostly enriched by large ellipticals. Poor clusters are indeed seen to harbor smaller 
amount of iron than rich clusters (e.g., Mulchaey et al.1993). 

In any case, it is almost certain that further X-ray observations of metal lines, espe- 
cially with ASCA, will be useful in elucidating the origin of the metallicity of the ICM. 

5. CONCLUSION 

We have discussed the hypothesis that gas ejected from dwarf galaxies in clusters 
can explain the observed intracluster gas. We have derived analytical expressions for the 
fraction of ejected gas, and the metallicity of the gas, and found good agreement with 
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previous numerical estimates. We have then calculated the metallicity of the resulting 
intracluster gas, and found that the hypothesis of dwarf galaxies as sources of the gas, is 
tenable only for clusters with very low metallicity. Additional sources, such as galactic 
winds from large ellipticals, are therefore needed to explain the existence of metals. We 
have not included the gas and metal mass from large ellipticals in our calculations, but 
only showed the deficiency of the hypothesis of dwarf galaxies as sources of intracluster 
gas, as this work is intended as a critique of this hypothesis. 

We further point out that, in CDM models, the gas in small-scale structures, which 
is dispersed and incorporated as diffuse intracluster gas by the process of hierarchical 
clustering, can account for most of the observed hot gas. We suggest that this process is 
of particular relevance in explaining the origin of the intracluster gas, especially for rich 
clusters. 



We are indebted to N. Arimoto, N. Trentham, and Y. Yoshii, for valuable discussions. 
BBN was supported by a fellowship from IUCAA. 
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Figure Captions 

Figure 1. (a) Fraction of gas mass 7 = M 9 /Mt expelled from galaxies as a function 
of the present-day galaxy mass (77 = 0.1, g = 0.1, r = 1.0). Thick solid curve: for an 
IMF with slope index x = 0.95. Thick dashed curve: for x = 1.35 (all curves are for 
Mi = 0.1, Mz )Sn = 8, M u = 100M Q ). (b) Iron abundance of gas Z Fe /Z &}Fe at the epoch of 
galactic wind as a function of the present-day galaxy mass (x = 0.95, g = 0.1, r = 1.0). 
Thick solid curve: for the fraction 77 = 0.1 of the supernova energy that is finally converted 
into the gas. Thick dashed curve: for r\ = 0.05. 

Figure 2. (a) Constraints on two parameters a (the slope of the dSph luminosity function 
at faint end) and f3 [(M/L) oc L 13 ] to reproduce the observed hot gas in Coma by gas 
ejection from dSph galaxies (/150 = 1, r\ = 0.1, g = 0.1, the upper mass bound of dSph is 
M + = 10 n M o ). Thick solid line: for the lower mass bound of dSph being M_ = 10 4 M Q 
and IMF index x of 0.95. Thick dashed line: for M_ = 10 4 M Q and x = 1.35. Thin solid 
line: for M_ = 1O 6 M and x = 0.95. Thin dotted line: when the uncooled primordial 
gas in the mass range Mu m = 1O 2 M ~ M_ = 1O 6 M given by eq.(17) is added to the 
expelled gas plotted as thin solid line, (b) Predicted iron abundance of ICM in Coma, 
Z^ e /Z QjFe , if all of the gas mass is enriched within dSph galaxies (x = 0.95, g = 0.1, 
M+ = 10 n M G ). Thick solid curve: for M_ = 1O 4 M and n = 0.1. Thin solid curve: for 
M_ = 1O 6 M and 77 = 0.1. Thin dashed curve: for M_ = 1O 6 M and 77 = 0.05. 

Figure 3. Fraction of uncooled gas erased from low- mass objects Mu m < M < M_ = 
1O 6 M in Coma (77 = 0.1, g = 0.1, x = 0.95, M + = 10 n M G ). 
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